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ABSTRACT 

Dwarf spheroidal galaxies are the most common type of galaxies, and are the most dark matter 
dominated objects in the Universe. Therefore, they are ideal laboratories to test any dark matter 
model. The Bose-Einstein condensate/scalar field dark matter model considers that the dark matter 
is composed by spinless-ultra-light particles which can be described by a scalar field. This model is 
an alternative to the A-cold dark matter model. In this work I study the kinematics of the dwarf 
spheroidal satellite galaxies of the Milky Way and Andromeda, under the scalar field/BEC dark 
matter paradigm in two limits: when the self interacting parameter is equal to zero, and when the self 
interacting parameter is 1. I find that dwarf spheroidal galaxies with very high mass-to-light ratios 
(higher than 100) are in better agreement with an NFW mass density profile. On the other hand, 
dwarf spheroidal galaxies with relatively low mass-to-light ratios and high luminosities are better 
described with the SFDM model. Such results are very encouraging to further test alternative dark 
matter models using the dynamics of dwarf galaxies as a tool. 

Subject headings: dark matter — galaxies: dwarf — galaxies: kinematics and dynamics — Local 
Group 


1. INTRODUCTION 

Dwarf spheroidal (dSph) galaxies are believed to be 
the most common type of galaxies in the Universe, and 
the building blocks of more massive galaxies in hierarchi¬ 
cal formation scenarios.The large mass-to-light ratios of 
these dSph galaxies suggest that the dynamics of such 
galaxies are dominated by dark matter (DM). There¬ 
fore, the study of dSphs is crucial for a better under¬ 
standing of the nature of DM. In particular, the dSph of 
the Local Group (LG) of galaxies satellites of the Milky 
Way (MW) and Andromeda galaxies, are good targets 
for study, since individual stars can be resolved and evo¬ 
lutionary histo ri es ca n be derived in great detail. 

1 Walker et all (l2QQ9h proposed that all dSphs follow a 
universal mass profile, which means that all dSphs are 
embedded within a universal dark matter halo. They 
found the best-fit for the relationship between velocity 
dispersion and half-light radius for a maximum velocity 
Vmax = 13 km/s, and a scale radius, ro = 150 pc, for 
a cored DM halo. For the NFW profile, the best fit is 
achieved for Vmax = 15 km/s, and a scale radius ro = 
795 pc. 

More recently, iCollins et all ([2014( 1 extended the pre¬ 
vious work, incorp oratin g kine matic information of 25 
Andromeda dSph (| Collins et al.ll2013( l. They found a dis¬ 
crepancy between the DM density profile best-fit for the 
dSphs in the MW compared with that for Andromeda’s 
dSphs (for both, cored and NFW cases). A good agree¬ 
ment between both MW dSphs and Andromeda dSphs 
samples is only reached, when three Andromeda outliers 
(And XIX, And XXI and And XXV) and three MW out¬ 
liers (Hercules, CVnI and Sagittarius) are removed from 
the dwarf galaxy sample. 

The discrepancy between the cuspy density profiles of 
DM halos predicted in simulations, with the cored den- 
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sity profiles derived from observations of dSph galaxies 
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H21 1 (the cusp/core 
problem), the overpopulation of dark substructure 
( Klypin et al. and the to-big-t o -fail problem 

( R,ead et al.l 2006c iBovlan et all I2011L I2Q12I ). have 
motivated alternative DM candidates to the ACDM 
model. 

An alternative model that lately has gained in¬ 
terest, is to consider that the DM is made of 
bosons described by a real (or complex) scalar 
field <f>: the Bose-Einstein cc mden s ate/scalar field 
DM model (BEC/SFDM) " 
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B EC/S FDM model is also known as Fuzzy DM 
(iHu et al.l l2000h or recently as Wave DM (jSchive et all 
120141 ). 

A very interesting feature of the BEC/SFDM model, 
is that it naturally produces cored halos. The size of the 
cores of such halos (for a fixed self-interacting parameter 
A) depends on the mass of the BEC/SFDM boson and 
the mass of the SFDM halo ( Mdm )• 

The SF DM model has been proved to be very 
successful (Mag ana fc Matosl 120121 ). It is consis¬ 
tent with the anisotropies of the cos mic mic rowave 
background radiation (CMB) (|Rodrfgnez-Montova et al.l 
12010( 1. and it can also achieve a better fit to high- 
resolntio n rotation curves of low-surface-brightness 
galaxies ( Robles fc Matos~ll2012[) , compared to the NFW 
([Navarro et al.l 1 1996 ) profil e. 

Lately, ILora et ahl ( 20121 ) and iLora fc Maganal ([2014( 1 
used the internal stellar structures of dSph galaxies to es¬ 
tablish a preferred range for the mass of the bosonic 
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particle. They performed TV-body simulations and ex¬ 
plored how the dissolution time-scale of the cold stel¬ 
lar clump in Ursa Minor (UMi) and Sextans depends on 
77 i 0 . They found that for a boson mass in the range of 
(3 < ? 7 i 0 < 8 ) x 10 -22 eV, the BEC/SFDM model would 
have large enough cores to explain the stellar substruc¬ 
ture in dSph galaxies. 

Moreover, iDiez-Teiedorl ( 2014 ) find a preferred scale 
radius of ~ 0.5 — 1 kpc, from the kinematics of the eight 
brightest dSphs satellites of the MW. 

In this work I investigate the idea of a univer¬ 
sal mass p rofi le fo r the dSph population of the MW 
(IWalker et al.l 2009 ) and Andromeda (M31) under the 
SFDM m odel , and compare it with a cored and NFW 
profile dCollins et al.ll2014t ). 

The article is organized as follows: in §2] the SFDM 
model is described and the Schrodinger-Poisson system 
is briefly reviewed. In 0 I discuss the universal NFW 
and cored DM profile for the dSphs in the MW and M31. 
In 0 I describe our results, and finally, in section §5] I 
discuss the results and give our conclusions. 


Qrl \3 2 . 

* = + U ^ + A l^! 2 ^ > ( 6 ) 

^(fU)=fU*- ( 7 ) 

The solutions of the Schrodinger-Poisson system are 
obtained assuming an harmonic behavior of the scalar 
field, such that 

t }{r,t) = e~( l ^(j>(f) , (8) 

where 7 is an adimensional frequency. 

Substituting 7/7 in Equations [6] and [7J the Schrodinger- 
Poisson system now reads as 


— M) = 2 f(C/- 7 ) + 2fA0 3 , 
dr z 



ref) 2 . 


(9) 

( 10 ) 


2. THE BEC/SFDM HALOS 

The DM halos can be interpreted as BEC/SFDM New¬ 
tonian gravitational configurations in equilibrium, which 
can be described by the so called Schrodinger-Poisson 
system, which is the Newtonian limit of the Einstein- 
Klein-Gordon equations 


,h M = 


H ' V 2 tp + Um«,tp+^\i:\ 2 <P , ( 1 ) 


2 m^, 


2 m^, 


W 2 U = 4t rGmJ#* 


( 2 ) 


The critical mass for such configurations is M cr u ~ 
0.6 ~ 10 12 Mq dR.uffini fc Bonazzolal 1969 : 

_ \ m f/> / _ 

iSeidel fc SuenlH991h . 

In Equations [T| and [2j m$ is the mass of the boson 
associated with the wave function 7 p. U is the gravita¬ 
tional potential (produced by the mass density source 
p = tti^It/I 2 ), and A is the self-interacting parameter. 
It has been shown that the la tter ter m dete rmine s the 
compactness of the structure dGuzman fc Urena-Lopezl 
2006 b 

It is convenient to work with the dimen sionless 
Schrodinger-Poisson system, given as follows ( Bernalll 
12001 


. d^ 1 y7 2 ? 


Uip + A|- 0| 2 -0 , 


( 3 ) 


V 2 U = ipip* , 


where 


A = 


A 


47rm^,G 


( 4 ) 

( 5 ) 


The solutions to the Schrodinger-Poisson system rele- 
vant for this work are those w ith spherical equilibrium 
(|Chi-Wai Laill2004tlBernalll2007D . Hence, it is convenient 
to work with the Schrodinger-Poisson system in spherical 
coordinates 


The SFDM halos are constructed by obtaining ground 
state (stabl e) solution s of the Schro dinger- Poisson sys¬ 
tem (I9lfl0l) (Guzman & Urena-Lopez 2004; B erna 3 Hqo 1 
iLora et ahll2012t iLora fc Magana 120141) . To guarantee 
regular solutions, the boundary conditions must satisfy 
that at r = 0 d r U = 0, d r <j) = 0, and that (f)(0) = <j) c =, 
where (f) c is an arbitrary value. 

The mass of this BEC/SFDM halo can be estimated 
as 

nOO 

M = 4tt (f) 2 r 2 dr , (11) 

Jo 

and must be finite number. The radius of this configura¬ 
tion is defined as 7 * 95 , which is the radius containing 95% 
of the mass. Note that both properties, the mass and 
the radius of the SF halo depend on the boson mass, and 
the self-interacting term. Three parameters </> c , and 
A, define a model completely. 

2.1. The A^> 1 case 

When the self interaction between the bosons in a 
Bose-Einstein condensate is taken into account, the 
Schrodinger equation can be interpreted as the mean- 
field approximation at zero temperature of the Gross- 
Pitaevskii equation (see Equation [Tj) . The limit where 
the number of particles in the BEC is very large, and thus 
the self interacting term dominates, is called the Thomas - 
Fermi limit (TFL) (|Pitaevskiill 196 It iDalfovo et al.lll999D . 

As the number of particles in the BEC be¬ 
come^ infinite, the TFL approximation becomes exact 
(jBarcelo et al.ll200~5l ) , giving as a result the classical limit 
of the theory. Then, the equations describing the static 
BEC in a gravit ational field with potential U take the 
following form dBohemer fc Harkolll2007 b 

vr(—) =i-pv(~) , ( 12 ) 

\m^J \mj 

V 2 U = AnGp . (13) 

With appropriate boundary conditions, one can in¬ 
tegrate the latter equations, assuming an equation of 
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Fig. 1.— Mass of the SFDM halo, as a function of the SFDM core radius for the A = 0 case. Each diagonal line corresponds to a 
different value of the SFDM boson m^\ 5 X 10 -23 (green), 10 _22 (blue), 2 x 10 _22 (pink), 5 X 10 _22 (red), and 10“ 21 (yellow) eV. The shaded 
region shows the perm itted core radius r c , such that the stellar clumps in Ursa Minor and Sextans are not destroyed (ILora et al.ll 12012: 
ILora & Maganal 120141b 


state of the form P = P(p). Assuming the non¬ 
linearity in the Gross-Pitaevskii equation of the form 
g{p) = ap T (where a and T a re constants greater than 
zero) (jBohemer fe Harkol 112007 ), then the BEC equation 
of state is given by 


P = P(p) = a(T - 1)/ . (14) 

If T is represented by T = 1 + ^ (where n is the poly¬ 
tropic index), then the equation of state of the static 
gravitationally bounded BEC is described by the Lane- 
Emden equation 


LA. 



(15) 


In particular, for a static BEC in the TFL with a 
politropic index n = 1, the Gross-Pitaevskii equation 
reduces to a Lane-Emden equation of the form 


d 2 0 2 dO n ^ 

U72 +737+^ = ° 

dtp f df 


with an analytical solution 

m = 


sin£ 

— 


(16) 


(17) 


Using the solution given by iBohemer fc Harko (2007 ) 
for the density profile, one has 


p(r) 


sin^nr / Rmax) 
P® nr/Rmax 

0 


if r < Rmax 
if T > Rmax i 


(18) 


where Rmax can b e interpreted as the size of the 
core radius of the SFDM halo. 


Finally, from Equation |TT] one can obtain the mass 
within a radius r in the TFL 


M(r) 


4 PpRmax 
7 r r 



irr 


Rr> 


-COS 



(19) 

Thus, the quantities that describe a system completely 
in the A 1 limit (TFL) are the central mass density 
p(c) and the size o f the S FDM halo R m ax- 
As iDiez-Teiedorl (2014D point out, the size of the DM 
halos in the TFL can be expressed as 


R 


max 


48.93 



2 


showing the dependency of R m ax with (yr %)- 


( 20 ) 


3. NFW AND CORED PROFILES FOR THE MILKY WAY 
AND ANDROMEDA DSPHS 

The measurements of central velocity dispersions 
of dSph galaxies very well constrain the dynamical 
mass within the deprojected half-light radii (M x / 2 ) 
(Wal ker et al. 2009; Wolf et aLll2010f h then one can com¬ 
pare such dSph mass (Mi/ 2 Vs r 1 / 2 ), with a certain in¬ 
tegrated DM mass profile, given a total DM halo mass 
M _ 

I Walker et all ( 2009 ) applied the Jeans equation to the 
velocity dispersion profiles of eight of the brightest dSph 
galaxies of the MW and consider the hypothesis that all 
dSphs follow a universal mass profile, i.e. that all dSphs 
are embedded within a universal dark matter halo. 

Using the density Equation, 

/ X __ Po _ 

P r (r/r 0 ) a (l + r/ro)P~ a 


(21) 
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IWalker et al.l (120091) fit a co red (a = 0 , /? = 3) and 
an NFW ([Navarro et al.l H996) (a = 1 , f3 = 3) DM pro¬ 
file. For the cored DM halo, they found the best-fit for 
the relationship between velocity dispersion and the half- 
light radius, for a maximum velocity V ma x = 13 km/s, 
and a scale radius, Rs = 150 pc. For the NFW profile, 
the best fit is achieved for V ma x = 15 km/s, and a scale 
radius Rs = 795 pc. 

ICollins et al ] (l20Ti extended the idea of a universal 
mass pro file, no w inc luding the M31 objects into the 
a nalysis (IColl ins et al.1 l2013h . 

ICollins et al.l (j2014f ) also use the NFW and core mass 
density profiles. They obtain the best fit for the NFW 
(cored) profile, to the whole population of dSphs, of 
Vmax = 14.7 =b 0.5 km/s (Vmax = 14.0 zb 0.4 km/s) and a 
scale radius of Rs = 876 zb 284 pc (Rs = 242 zb 124 pc). 
They concluded that neither value is a good fit for many 
of the LG dSphs. They found that if they do not 
take into account five outliers, three in M31 (And XIX, 
And XXI, and And XXV), and two in the MW (Her¬ 
cules and CVn I), the fit is significantly better for the 
NFW (core) mass profile, V max = 16.2±f;f km/s (V rnax = 
15.6±iJ km/s) and R s = 664t^| pc (R s = 225+ 5 ° pc). 

4. A UNIVERSAL SFDM HALO? 

4.1. The A = 0 case 

One can consider the self-interacting term A to be 
zero in Equation [T| (i. e. the self-interaction is negli¬ 
gib le), such a ca se is also known as the fuzzy DM model 
( Hu et ahl 1200(f) . 

In Figure [TJ the mass of the SFDM halo is plotted as 
a function of the core radius (r c ) for A = 0. The core 
radius is defined as the radius at which the initial density 
has dropped by a factor of two. 

When analyzing the behavior of the different 
models with fixed Mdm and ra^, I found out that 
there is a correlation between the mass of the 
SFDM halo and the core radius, for a given m^. 
The mass of the SFDM halo is inversely propor¬ 
tional to the DM core radius independently of 
the value of (i. e. Mdm ck AT 1 )? but the inter¬ 
cept varies depending on the value of m^. Then, 
I found a fit over different values of and their 
corresponding intercepts (o), and found the rela¬ 
tion between the intercept and the value of of 
a oc m 8,T . 

The complete relation between the mass of the dark 
matter halo and the dark matter halo’s core radius for a 
given m^, is given by 

M dm = lO^r- 1 , (22) 

where a depends on the mass of the boson as 

= 10“ 13 * 9 ar 8 ' 7 . (23) 

For a given Mdm and m^, there is a slope o, for which a 
unique r c will be defined. With the latter relations, one 
can build a set of continuous models for a given Mdm 
and m^. The continuous models are shown as diagonal 
colored lines in Figure [lj Each color represent a mass 
of the boson (green 4 5 x 10 -23 eV, blue 4 1 x 
10 -22 eV, pink 4 2 x 10 -22 eV, red 4 5 x 10 -22 eV, 
yellow 4 1 x 10 -21 eV). 


It is clear to see that for a given halo mass, lower values 
of m/s have larger core-radius, and higher values of mn$ 
have smaller core-radius. For example, for a fixed SFDM 
halo mass of 10 8 M©, if = 10 _22 eV the correspond¬ 
ing core radius is ^ 3.5 kpc. But if m$ = 10 -21 eV the 
corresponding core radius is ~ 30 pc. 

I consider three different groups of dSph data; the 
dSphs in the MW (19 galaxies), the dSph in the M31 
(22 galaxies), and al l the galaxie s (41 galaxies). In 
order to compare with ICollins et al.l ( 2014 ) results, I do 
not take into account five dSphs (Sagittarius, Hercules, 
CvnI, AndXIX, AndXXI and And XXV). 

In order to see how well these three groups of dSph 
galaxies can be fit with a single SFDM profile (for a 
SFDM halo mass M and a SFDM boson mass m/, I 
use the following maximum likelihood fitting routine 


LsFDM({rh,i, V C: i,8y c ^}\m ( f ) ,M) = 


n 



x exp 


(Vc, 


SFDM 


- v c ,y 




(24) 


where V c ,sfdm is the circular velocity predicted by the 
SFDM model, rh,i is the half-light radius of the i— dSph, 
V c s is the measured circular velocity at the half-light 
radius, and 6y c ,j \s its uncertainty (the observationa l data 
are t aken from iMcConnachiel I 2012 L ICollins et al.l 120141 
and lMartin et al.ll2015l ). 

Equation [24] is a measure of how well a model with 
specific parameters and M represents a given group 
of galaxies. In order to find the best fit for the three 
samples I use the Evolution Strategy with Covariance 
Matrix Adaptation (CMA-ES) optimization. For details 
on the CMA-ES routine, see Appendix [A] 

The best fit for the 19 MW sample is M = 5.12 x 
10 7 M© and a = 3.8 x 10 -22 eV. For the 22 M31 
sample the best fit corresponds to M = 4.44 x 10 7 M© 
and a = 4.53 x 10 -22 eV. For the whole sample (41) 
the best fit corresponds to M = 4.8 x 10 7 M© and = 
4.17 x 10 -22 eV (see Tabled]). The MW fit gives as a 
result a slightly more massive DM halo, and a smaller 
mass of the boson, compared with the M31 and the all 
data sets. As a consequence the core radius for the MW 
sample, is the largest (r c ~ 0.45 kpc). Such a size of the 
DM core radius is large enough to ex plain t he l ongev ity of 
the old-cold stellar clump in UMi (iLora et al.ll2012[), and 
the tw o stellar substructures in Sextans (jLora fe Maganal 

Uni). 

In the upper panels of Figure [2j I show the half-light 
radius circular velocity (V Cj i /2 = ) as a func¬ 

tion of the half-light radius of all M31 and MW dSphs 
in the sample, along with the best fit to all sample 
for an NFW density mass profile (see dashed dark-blue 
line) an d a cored de nsity m ass profile (see dashed purple 
lines) (ICollins et al.l !'2014lh and the SFDM profile (see 
light-blue dash line). The shaded regions in Figure [ 2 ] 
indicate the la error. The computed maximum circu¬ 
lar velocity for the SFDM model for the all sample is 
Vmax ~ 16.63 km/s. This value of the maximum circular 
velocity is in a very good agreement with both computed 
values of the NFW profile (Vmax = 16.2+ 2 ;y km/s), and 
the core profile (V ma x = 15.61/J km/s). 
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TABLE 1 

Parameters of the best fit after a CMA optimization to the SFDM profile for three groups: the MW dSph(19), the M31 

dSph(22), and ALL the sample(41). 


Group 

Number 
of galaxies 

M 

[10 7 M 0 ] 

m c/> 

[10- 22 eV] 

e 

icr 5 

^95 

[kpc] 

V c 

[km / s] 

r core 

[kpc] 

M31 

22 

5.12 

3.80 

7.07 

0.90 

16.15 

0.46 

MW 

19 

4.44 

4.54 

7.31 

0.73 

16.71 

0.37 

ALL 

41 

4.80 

4.17 

7.28 

0.80 

16.63 

0.41 



h /2 [pc] 


• M31 






h /2 [pc] 



Fig. 2. — Best fit to the NFW (dark blue), and the core (purple) DM density profile from [Collins et ahl (120140 . leaving out the MW 
dSphs’ outliers (Hercules, CVnI and Sagittarius, see white triangles), and three M31 outliers M31 (AndXIX, AndXXI, AndXXV, white 
circles). The SFDM profile studied in this work, is shown in light-blue. The shaded regions correspond to a lcr deviation, from each DM 
model fit. The upper panels show the circular velocity as a function of the half-light radius, for ALL, M31, and MW samples, respectively. 
The lower panels show the DM halo mass, as a function of the half-light radius, for ALL, M31, and MW samples, respectively. The dSph 
galaxies used in this work are over plotted with their corresponding uncertainties. The gray circles correspond to the M31 dwarfs, and the 
black triangles correspond to the MW dSphs. 


The dSph galaxies with low V c s and small half-light 
radii can be well reproduced with the cored and the NFW 
density profile, but cannot be well reproduced with the 
SFDM profile. On the other hand, the dSphs with V c s 
between 6 — 10 km/s, and high values of the half-light 
radii 1000 pc) can be reproduced with the SFDM 
model, but cannot be well reproduced with either the 
NFW profile, or with the core profile. 

In the lower panels of Figure [2j I show the same fits 
(core, NFW and SFDM) for the mass of the DM halo, 
as a function of the half-light radius. It is very clear 
from this Figure, that the NFW fits the all data better 
at small-mass-small-half-light-radius, whereas the SFDM 
fits better the all data at high-mass-high-half-light- 


radius, where both core and NFW profiles fail. 

In order to look in more detail at the properties of the 
dSphs that are in better agreement with the SFDM pro¬ 
file, I plot the mass inside the half-light radius as a func¬ 
tion of the luminosity at the same radius (see left panel 
of Figured]). The diagonal dashed lines in Figure [3] cor¬ 
respond to four fixed mass-to-light ratios ( M/L = 1000, 
100, 15 and 3). The circles represent the M31 dSphs, and 
the triangles represent the MW dSphs. 

For each dSph I computed the ratio £ between the theo¬ 
retical mass values from the best DM fi ts of the all dat a 
set, and the observed mass values ([Collins et al.l 12014). 
( = Mdm /MdSph- I then define that a dSph is well re¬ 
produced by a certain DM model if ( is at most 2. The 
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latter means that the dSph mass obtained from the best 
DM model differs from the observed value by a maxi¬ 
mum of a factor 2. A dSph galaxy is in agreement with 
an specific DM model when ( is minimized. 

In the left panel of Figure [3j the color code is as follows: 
the light-blue symbols (both M31 circles and MW trian¬ 
gles) represent the dSphs which are in good agreement 
with the SFDM profile. The purple symbols show the 
dSphs which are better reproduced with the cored DM 
profile, and the dark-blue symbols represent the dSphs 
which are in better agreement with the NFW profile. On 
the other hand, the yellow symbols correspond to the 
dSph that are well reproduced with the three DM mod¬ 
els, whereas the red symbols show the dwarf galaxies that 
cannot be reproduced with any DM model studied in this 
work (i. e. £ > 2). 

One can immediately see from Figure [3] that the dSphs 
which are in a better agreement with an NFW DM pro¬ 
file are very dark matter dominated, with mass-to-light 
ratios between 1600 > M/L > 100. There are only two 
exceptions Leol and LeoII, which mass-to-light ratios of 
7.1 and 12.9, respectively. 

The galaxies which are better reproduced with an 
NFW profile are also the less luminous of the sample 
with L < 10 4 L 0 . The only outlier is UMal, which is 
in a better agreement with the SFDM model. On the 
other hand, it seems that these galaxies have a mass at 
the half-light ratio limit of ~ 10 7 M 0 . 

The upper-right gray shaded region in Figure [3] con¬ 
tains a set of five M31 dSphs, and two MW dSphs which 
are in good agreement with the SFDM profile (light-blue 
symbols). If one also takes into account the yellow sym¬ 
bols, which are dSphs well reproduced with either DM 
model, then four more can be added, resulting in a set 
of 11 galaxies which are very well reproduced with the 
SFDM model. These galaxies share the property of being 
less DM dominated than those which are better repro¬ 
duced with an NFW DM profile. The mass-luminosity 
ratios of these galaxies range from ~ 100 to 5. On the 
other hand, these dSphs are the most luminous galaxies 
of the sample, with a luminosity ranging from ^ 10 5 to 
5 x 10 6 L 0 . It is remarkable that there are no dSph in 
this region which are more in agreement with an NFW 
profile. 


4.2. The A >> 1 case 

In this subsection, the kinematics of the dSphs in the 
MW and M31, for large values of the self-interacting pa¬ 
rameter A , are analyzed. In the A 1 limit, as men¬ 
tioned in section P~2l the SFDM density profile is de¬ 
scribed by Equation [18] In order to build the DM halos 
the value of the central DM density po and the size of 
the SFDM halo Rmax must be fixed. 

To study how well the circular velocity, and mass pro¬ 
file of the three groups of dSphs (M31, MW and ALL) 
are fit with a single SFDM halo in A >> 1 regime the 
maximum likelihood fitting routine is modified, so that 
po and Rmax are the free parameters. Then similarly to 
Equation [24] one has 


L T F({rh,i, V Cti , Sv c ,i}\po, R max ) — 


n 



x exp 


{Vc,TF - 

2< ^Vc,i 


(25) 


where W,tf is the circular velocity in the TFL. 

Using the CMA-ES optimization routine described in 
Appendix [3 I find that for the M31 dwarfs, the best fit is 
achieved for a value of m ^/A = 3180, which corresponds 

to a po ~ 8 x 10 7 M 0 kpc -3 , a Rmax — 0.87 kpc, and 
at circular velocity value of 18.24 km/s (see Table [2]). 
The latter value is approximately 13% higher than the 
one obtained for the A = 0 case (16.1 km/s). The core 
radius (halo DM mass M) in the TFL for M31 dwarfs, 
is 0.52 kpc (6 x 10 7 M 0 ), very similar to the 0.45 kpc 
(5 x 10 7 M 0 ) find for the A = 0 case. The results are 
summarized in Table [2} 

For the MW dwarfs the central density is the highest of 
the three groups, 1.8 x 10 8 M 0 kpc -3 , whereas the maxi¬ 
mum DM radius for the MW dSphs is the smallest of the 
three groups, Rmax = 0.59 kpc (which corresponds to a 
core radius of 0.36 kpc). Such a core radius is very sim¬ 
ilar to the one obtained for the A = 0 case, and thus, it 
is large enough to explain the stellar substructures found 
within the UMi and the Sextans dSphs. The maximum 
circular velocity is 18.6 km/s, again approximately 11% 
higher than that found in the A = 0 case (16.7 km/s). 
For the MW dwarfs a value of m ^/A = 6788 is obtained 
(i. e. m^/ A 1 / 4 « 9). 

The upper panels of Figure [4] show the circular velocity 
as a function of the half-light radius, and DM halo mass 
as a function of the half-light radius (lower panels) for 
the SFDM model for A = 0 and A >> 1. The pink 
lines show the best fit for the SFDM model in the TFL 
(A >> 1) for each of the dSph data sets: ALL (left 
panels), M31 (central panels) and MW (right panels), 
respectively. The shaded regions show the la deviation. 
The results obtained for the A = 0 case (light-blue), are 
also plotted in order to facilitate comparison between 
both SFDM profiles. 

In the right panel of Figure [3] I show the seven dSphs 
which are in better agreement with a SFDM model in the 
TFL (such galaxies are plotted with pink symbols) also 
including the previous studied DM profiles. The galaxies 
that are better reproduced with the SFDM in the TFL 
are also mostly contained in the upper-right shaded gray 
region. This means that the SFDM model in the TFL 
better fit galaxies with low mass-to-light ratios. 


4.3. The universal SFDM profile compared 

In lLora et al.l (2012 ) we performed 7V-body simulations 
of the UMi dSph,^which contains a cold-old stellar sub¬ 
structure (jKlevna et al.ll2003f ). We explored how the dis¬ 
solution time-scale of such stellar substructure depends 
on the mass of the boson ra^. The boson mass range 
obtained from UMi’s dynamics is 0.3 x 10 -22 < < 

10 -22 eV. 

Similarly, in iLora fc Maganal (2014) we investigated 
the dSph Sext ans, which has two different old-cold stella r 
substructures ([Battaglia et al.ll2011l : I Walker et al.ll2006[ ). 
For Sextans we require a mass for the SFDM boson of 
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Fig. 3.— Mass at the half-light radius as a function of the luminosity at the half-light radius. The M31 dSphs are shown with circles 
and the MW dSphs are shown with triangles. The dark-blue(purple) symbols show the dSphs that are in better agreement with an NFW 
(cored) DM profile. The light-blue symbols show the dSphs which are in better agreement with the SFDM profile. The yellow symbols 
show the dwarfs that are well reproduced with the three DM models, and the red symbols show the dwarfs that are not well reproduced 
with any DM model. The white symbols represent the galaxies that were not taken into account (Hercules, CVnI, AndXIX, AndXXI and 
AndXXV). The diagonal dashed lines show the mass-to-light ratios M/L = 1000, 100, 15, and 3. In the right panel, the dSph that are in 
a better agreement with the SFDM in the TFL i. e. A >> 1, are also plotted (see pink symbols). 


TABLE 2 

Parameters of the best fit in the TFL (A >> 0) for three groups: the MW dSph(19), the M31 dSph(22), and ALL the 

sample(41). 


Group 

Number 
of galaxies 

po 

[10 7 M 0 kpc 3 ] 

Rmax 

[kpc] 

Mrnax 

[10 7 M@] 

[eV 4 ] 

AV4 

[eV] 

Vc 

[km/s] 

r core 

[kpc] 

M31 

22 

8.06 

0.86 

6.716 

3180.96 

7.51 

18.24 

0.52 

MW 

19 

17.96 

0.59 

4.80 

6788.80 

9.08 

18.64 

0.36 

ALL 

41 

9.74 

0.85 

7.60 

3325.20 

7.59 

19.61 

0.51 


0.12 x 10 -22 < rrify < 8 x 10 -22 eV, in order to guarantee 
the survival of both stellar substructures. 

In Figure [TJ the corresponding mass of the boson is 
plotted (purple blue stars) for UMi and for Sextans. The 
error-bars denote the la deviation. The best fit of the 
MW(triangle), M31 (circle), and all the sample (square) 
are also plotted. 

For the Sextans dSph, all the values for the mass M 
and the core radius r c in the light-gray area in Figured] 
are permitted. In particular, those corresponding to the 
three best SFDM fits (MW, M31 and all ). The lat¬ 
ter means that the dynamics of the Sextans dwarf, is in 
agreement with a unique SFDM mass profile. The UMi 
dSph is located above the MW, and all fits in Figured] 
being the most MW dSph restrictive case. However, it 
is still in agreement with the MW (circle) case at la 
confidence. 

My findings in this work in the TFL can be directly 
compared with the recent results of lUiez-Teiedorldml . 
They constrained the parameters of the SFDM self¬ 
interacting parameter with the kinematics of the eight 
brightest dSph of the MW. They reported a preferred 
value of Rmax ~ 1 kpc, which corresponds to m$j A 1 / 4 ~ 


7 eV. In this work I obtain for the MW sample a value of 
Rmax ~ 0.6 kpc, which co rresp on ds to m^/A 1 / 4 ~ 9 eV. 
It has to be noted that IDiez-Teiedorl ( 2014 ) analyzed 
the eight brightest dSph in the MW; in this work the 
MW sample was comprised of 19 dwarf galaxies. If the 
MW sample is restricted to the eight brightest dSph, I 
obtain a value Rmax = 0.57 kpc, which corresponds to 
ra^/A 1 / 4 = 9.24 eV. This result is still in agreement with 
the results obtained from the 19 MW dwarf sample, but 
ra^/A 1 / 4 (and Rmax) is somewhat higher (lower) that 
t he mVA 1 / 4 ~ 7 eV (Rmax ~ 1 kpc) values reported by 
IDiez-Teiedorl (|2014t ). 

It has to be noted that iDiez-Tejedorl (j2014h an¬ 
alyze the velocity dispersion of the eight bright¬ 
est dSph galaxies in the MW, arguing that the 
stellar component of each of the eight galaxies 
are in dynamical equilibrium, and that the stel¬ 
lar distribution traces the DM distribution. They 
fit three free parameters in their models: Rmax, 
M ma x and the orbital anisotropy of the stellar 
component of each galaxy, which is a difference 
in the approach of this work. It is encouraging 
that, even with these differences, when the uncer- 
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Fig. 4.— Same as Figure [2] but for the best fits of the SFDM model, and the SFDM in the TFL. 


tainty averaged over the eight dSphs is taken into 
account in Rmax (and thus m^/A 1 / 4 ) one obtains 

Rmax = 1-05tg 22 (m <j> /\ 1/4 = 6.82t(J;^), which is 
roughly in agreement with the findings reported 
in this work. 

IStrigari et al.1 ( 2008 ) suggested that dwarf galaxies 
satellites of the MW, have a common mass within 300 pc 
of ^ 10 7 M 0 . They suggest that such a finding could 
shed light in a characteristic scale for the clustering of 
dark matter . In Jdgure [5] the gray shaded region shows 
IStrigari et al.l (2008)’s results. Overplotted are the re¬ 
sults of this work. The light blue symbols correspond 
to the A = 0 case, and the pink ones correspond to 
the A 1 case. It is encouraging to find that for the 
best fits of the SFDM model for the M31 and MW dSph 
galaxies, the mass within 300 pc is in agreement with 
IStrigari et al.l (|2008l ) findings. 

It has to be noted that the central mass density of 
th e MW’s dS phs r ange from po_~ 0-03 to 0.3 M 0 pc -3 
([Kormendv fc Freemanl 12014 iBurkertl I2014D . It is en¬ 
couraging that the latter values of the central density 
are not only in agreement with the central density ob¬ 
tained for the MW’s sample 0.18 M 0 kpc -3 ), but 
also with the whole sample 0.1 M 0 kpc -3 ). 

5. CONCLUSIONS 

The very high dark-to-stellar mass ratios of dSph 
galaxies suggest that they are the most DM dominated 
objects in the universe, and therefore ideal laboratories 
to test any DM alternative model. In this work I compare 
the kinematics of 22 dSph galaxies satellites of M31, and 
19 dSph galaxies of the MW. I study also a third group 
of galaxies containing the sum of all M31 and MW dSphs 
(ALL). 
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Fig. 5. — The gray shaded region shows the results of 
IStrigari et al.l ( 20081 "). The resulting mass within 300 pc for M31 
(circle), the MW (triangle), and all the sample (square) are shown, 
for the A = 0 case (light blue symbols) and the A 1 case (pink 
symbols). 

I study the hypothesis that all dSph are e mbedded 
in dark matter halos w ith a same mass (Wa lker et al.l 
120091: ICollins et al1l2014lh by fitting the dSph’s half-light 
radius, and the velocity at the half-light radius, to the 
SFDM model, in two different regimes: when A = 0 and 
when A 1. 

There is a very good agreement between the velocity at 
the half-light radius for the best fit for the ALL sample 
(16.4 km/s) from the A = 0 case, with those obtained 
using an NFW (16.2 km/s), and a cored (15.6 km/s) 
DM profile (see Equation ED). This corresponds to a 
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mass of the SFDM boson of ^ 4 x 10~ 22 eV, whic h is in 
agreement with our p revious findings (jLora et al.ll2Q12t 
iLora fc Maganall2014[) . 

A higher value for the velocity at the half-light radius 
was found for the A >> 1 case (19.6 km/s) resulting in 
a niff,/ A 1 / 4 ~ 7.6 eV. 

The M31-dSph galaxies with high luminosities and 
mass-to-light ratios ranging from 27 to 78, are better 
reproduced with the SFDM model than with an NFW 
or core DM model (see light-blue symbols in Figure [3]). 
Hence, MW-dSph with very high mass-to-light ratios 
(ranging from 1500 to ^ 100) and very low values of 
the luminosity are in better agreement with a NFW/core 
DM profile than with a SFDM (see dark-blue and purple 
symbols in Figure [3]). 

The mass within 300 pc for the SFDM model (for both 
A = 0 and 1) is in a very good agreement with 


iStrigari et al.l ( 200 8)'s findings. The latter suggests a 
universal SFDM halo mass for the dSph galaxies in the 
MW and M31. These results are encouraging, since also 
a unique mass of the SFDM boson (for the A = 0 case) 
of ^ 4.8 x 10 -22 eV is obtained. One would expect a 
unique mass of the SFDM boson, and not to have a dif¬ 
ferent SFDM boson mass for different dSph galaxies. 

I would like to thank Michelle L. M. Collins and Nico¬ 
las F. Martin, for making their data available to me. I 
also thank Steffen Brinkmann and Cohn W. Glass from 
HLRS, for making their CM A routine available to me. 
I thank Andreas Just, Juan Magana and Avon Huxor 
for very helpful comments and discussions, that gave as 
a result an improved version of this paper. I gratefully 
acknowledges support from the Heidelberg University In¬ 
novation Fund FRONTIER. 
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APPENDIX 

THE EVOLUTION STRATEGY WITH COVARIANCE MATRIX ADAPTATION OPTIMIZATION 

The Evolution Strate gy with Covariance Matrix Adaptation (CMA- ES) is a stochastic, derivative-free method for 
numerical optimization (jHansen fc OstermeierlDOOH lHansen et al.ll2003[) . It is robust and superior to derivative based 
search methods in rugged search landscapes, i.e. containing noise, outliers, large derivatives etc. 
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First Generation Second Generation Third Generation 


Fig. 6.— Representation of the CMA optimization method. 

An evolutionary strategy algorithm uses mutation and selection on a given sample of solutions until a termination 
criterium is met. More precisely, the algorithm is initialised with a population of (i candidate solutions in the solution 
space, called parents. From these parents, A > /i new candidate solutions (children) are created by mutation, i.e. by 
scattering their properties (=genes) around the parents properties: 


~ aJ\fi(m,C) for i = 1 ,..., A . (Al) 

In the latter Equation is the property vector of child i, m is the weighted mean property vector of the parents, a 
is the step size, A f is the normal distribution around m, which is distorted by the covariance matrix C. The children 
are evaluated by a fitness function (/(x) : M n —* M) and only the fittest /i children become the parents of the next 
evolutionary step. 

CMA-ES enhances the evolutionary approach by updating the step size and the covariance matrix independently 
after each evolutionary step. For that purpose, the algorithm takes into account the /i fittest children according to a 
given fitness function, the distance and direction of the new weighted mean relative to the old one, and lastly, how 
many children achieved a better fitness than the best parent. 

The mean of the distribution is updated such that the likelihood of previously successful candidate solutions is 
maximized. The covariance matrix of the distribution is incrementally updated such that the likelihood of previously 
successful search steps (mutations) is increased. Both updates can be interpreted as a natural gradient descent. See 
Figure [6J 

The adaptive algorithm design makes CMA-ES a black-box optimization strategy in the sense that no information 
about the solution landscape is required. The algorithm is furthermore robust against transformations of the solution 
space, e.g. by exchanging the order of parameters in the fitness function or different dimensional scaling of the solution 
space. Dependences between the solution candidate properties influence the covariance matrix, thus accelerating the 
search for the optimal solution. 





